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The topology, the symmetry involving the shape of dielectric cylinders, and the lattice structure 
are among the most important ingredients in the architecture of photonic crystals. In this paper, 
we present a systematic derivation of the formulas which are needed in computing the photonic 
band structures of many commonly used two dimensional lattice structures and dielectric cylinders 
with various kinds of symmetries and rotations. Further the results are applied to arrays of hollow 
cross-shaped cylinders embedded in the alumina ceramic background. A large complete photonic 
band gap is found in the high frequency regime. 
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PACS numbers: 41.10.Hv, 71.25.Cx, 43.20.+g. 



I. INTRODUCTION 

> 

oc 

' Over the past ten years, the propagation of classical waves in a periodic medium has attracted considerable interest. 
This includes electromagnetic (EM) wave propagation in periodic dielectric structures, and acoustic and elastic wave 
propagation in periodic elastic composites. A new research field emerges as the wave crystals including both photonic 
and sonic crystals. The photonic or sonic crystals respectively refer to crystal-like structures that modulate EM or 
acoustic wave propagation and thus lead to dispersion bands, in analogy with the electronic energy bands in solid 
^ . state physics [Q. 

The research on photonic crystals has been particularly intensified, after the suggestion that photonic band gaps 
(PBG) could hinder spontaneous emission and block propagation of EM waves, thus providing the possibility to 
manipulate the propagation of EM waves . Photonic crystals offer an unparallel opportunity to design new optical 
devices and hold a great potential for many significant applications, such as semiconductor lasers and solar cells [^[Q, 
high quality resonator and filters Q, controlling photon emission H, optical fibers [Q, guiding and bending of EM 
, waves with minimum losses |g-{[0), single mode waveguides for light [ Jllj , low dimensional efficient transport of electrons 
and excitons by nanostructural networks , all polymer optoelectronic devices jl3j , semiconductor memory cells [ pd| . 
Many methods have been proposed for fabricating photonic crystals. These include square spiral microfabrication 
achitecture for large three dimensional band gaps |15| , filling the voids in titania with air by precipitation for the optical 
spectrum [|l6f , using three dimensional carbon structures Jl^l , large scale synthesis of silicon photonic crystals [fl8| , 
fabrication of photonic crystals for visible spectrum by holographic lithography []l9|| , the electrochemical techniques 

IB 

Indeed, the past a few years have witnessed rapid advances in both better understanding of the exquisite properties 
of photonic crystals and manipulation of EM waves by photonic crystals. A rich body of literature on photonic crystals 
exists and can be found on the internet [^2| . Recently, scientists also investigated the spine from sea mouse [^3| . They 
discovered that the spine consists of an array of regularly arranged hollow cylinders, and this simple structure gives 
rise to a spectacular iridescence. This is a remarkable example of photonic crystals by a living oraganism. 

Although three dimensional (3D) photonic crystals suggest the most intriguing ideas for novel applications, two 
dimensional (2D) structures also find several unique uses (24), including the aforementioned waveguides and com- 
munication fibers @,|ll| and the 2D periodic structures in living animals p3[ , feedback mirror in laser diodes |25| ] 
and so on. In addition, fabricating 3D periodic structures in the near infrared regime poses a significant challenge 
compared to the two dimensional situations Jl7| , p6 -28|. Due to these reasons, the study of 2D photonic crystals has 



been overwhelmed in the last few years. The earliest theoretical analysis of 2D photonic band structure was done 
by Plihal et al. |2{| . The experimental observation of photonic band structure in 2D periodic dielectric arrays was 
subsequently reported by Robertson et al. ]3C| ]. 

The important issue in the fabrication of photonic crystals is to create large, robust complete band gaps within 
which propagation of EM waves is prohibited in any direction. Several methods have been suggested for obtaining 
large complete band gaps in 2D situations. For example, it has been shown that large band gaps can be obtained 
by such as varying dielectric contrast ratio and filling factors, inserting a third component into the existing photonic 
crystals |§jj ], reducing the structural symmetry |^| , using non-circular rods |}2 33 and subsequently by rotating 



the non-circular rods 34 3q|, rotating the lattice structures 36 1, using anisotropic dielectric materials ]37|] , using the 
effects of magnetic permeability p8| , using metallic or metallodielectric rods ]39|-|4l|l, placing rods of various shapes 
on different lattice configurations such as square |||] , triangular 37 4^] , honeycomb and so on. Each approach 
may have its advantages and shortcomings. For example, the dielectric contrast is limited by material availability. 
The symmetry reduction and using non-circular rods may reduce the degeneracy of photonic bands at high symmetry 
points in the Brillouin zone, thus increasing band gaps. Although the metallic photonic crystals can yield large band 
gaps, they suffer from absorption. While the symmetry reduction method can enhance some high order band gaps, 
the low order gaps are often reduced [[uj. Therefore each method has its own applicable situations. 

Inspecting these progresses made towards better design of two dimensional photonic crystals, we realize that the 
topology, the symmetry involving the shape of dielectric cylinders, and the lattice structure are among the most 
important ingredients in the architecture of photonic crystals. The different combinations of these factors lead to 
applications for various purposes. We are therefore led to the task of deriving systematically necessary formulas 
computing band structures for various configurations. This paper is one of our attempts. In this paper, in an 
organized fashion we present the analytic results for computing the band structures of most commonly used lattice 
structures and dielectric cylinders with many kinds of symmetries and rotations. With the aim in mind that the 
reader can readily make use of these results, we summarize them in tables. 

This paper is organized as follows. The general theory for EM waves in an arbitrary 2D periodic structures 
is presented in the next section. Using the plane wave expansion method, the secular equations are derived for 
determining the band structures for both E-polarization mode and H-polarization mode with the electric field and 
the magnetic field parallel to the longitudinal axis respectively. In the formulation, a structure factor is identified. 
Once this factor is known, the band structure can be computed by a standard diagonalization method. In section III, 
the structure factor is derived for a general configuration. Some common crystal structures with several rod shapes 
are considered in section IV. The results are presented in two tables. A numerical example is shown in section V, 
followed by concluding remarks in the last section where the extension to sonic crystals is also discussed. 

II. THEORY 

The photonic band structures can be obtained by solving Maxiwell's equations using the plane-wave expansion 
method. This method can be referred to, for example, Refs. [ p9||34|j4^ , ^| . For the sake of convenience, we provide a 
brief account of this method. Consider a periodic array of dielectric cylinders. The longitudinal axes of the cylinders 
are along the z-axis. For dielectric materials, Maxiwell's equations are represented in terms of the magnetic field H 
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-H(r), (1) 



where e is the position-dependent dielectric constant, r is the coordinates in the plane perpendicular to the rods, and 
Co is the EM wave phase speed in vacuum. By Bloch theorem jjj, the magnetic field can be written as 

H(r) = e lkr H fc (r), (2) 

in which H/^r) is a periodic function of the lattice structure, and k is the Bloch wave vector within the first Brillouin 
zone. 

Using the Fourier transformation, the magnetic field and the dielectric function can be expanded as, 

H *( r ) = E E e^,k(G)e iG -, (3) 

G 3=1,2 

and 

£ (r)=5>(GK G ' r , (4) 

G 

where ij is the base vector for the magnetic field, G is a reciprocal lattice vector. Hereafter, all quantities with a hat 
refer to the unit vectors. For two dimensions, the E-polarization and H-polarization modes are decoupled. By taking 
Eqs. (||), (||), and (|J) into Eq. (|l|), we obtain two eigen-equations 

V |k + G| |k + GV^G - G')ffx,k(G') = ^fli jk (G), (5) 
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for the E-polarization and 



£(k + G) • (k + G')e-\G - G')if||, k (G') = ^H ]]tk (G), (6) 

G' C 

for the H-polarization. Here e~ x (G — G') is the inverse matrix of e(G — G'), _L and || refer to the direction perpen- 
dicular and parallel to the z-axis. The Fourier components e(G — G') is calculated as 

e (G) = i^e(r) e - Gr dr, (7) 

where the integration is performed over the area of one lattice unit cell, and A is the area. 
For binary situations, Eq. (|7]) can be further simplified into 

e (r\-l + forG = 0, 

[ ' ~\(e a - e b )S(G) for G jt 0, ( *> 

with e a and e& referring to the dielectric constants for the cylinders and background separately, and / is the filling 
factor defined as fraction of the area occupied by the cylinders in one unit cell. The factor S(G) relies only one the 
geometry of the cylinders and the lattice structures, and is given by 

S(G) = - I e-* G 'dr, (9) 

A JA d 

where the integration is carried over the area occupied by the cylinders in the unit cell. We name S(G) as the 
structure factor. 



III. EVALUATION OF STRUCTURE FACTORS 



It is clear from the above derivation that once the structure factor is known, the band structure for either the 
E-polarization or H-polarization modes can be readily evaluated by taking the structure factor into Eq. (||) and 
subsequently into Eqs. (^) and ([|) respectively. In the following we evaluate the structure factor for common crystal 
structures. 



A. Transformation of S(G) under some operations 



To proceed, first we discuss how the structure factor changes under certain operations with regard to a single 
cylinder. Different cylinders in a unit cell may undergo different operations. The final structure factor of the unit 
cell will be the sum of each individual structure factor. These operations include translation, rotation, reflection and 
scaling (dilation and contraction). Under these operations the new 5(G) is related to the original one by applying 
simple transformations. The knowledge of these transformations are useful in the calculation of the structure factor 
for various lattices. 



1. Translation 



In the x — y plane, the translational operation is the simplest operation and it transforms the coordinates by a 
constant displacement r as 

T t [r] = r' =r + r . 
Under this operation, the structure factor changes as 

T t [5(G)] = S' = e- lG - r "S(G). (10) 



2. Rotation 



This operation rotates the cylinders by an angle 9. Since the calculation ol structure factor must be coordinate 
independent, we can rotate the coordinate system by the same angle 9 to simplify the calculation. In the new 
coordinates, the structure factor will be the same form as the original one. However, in the new system the vector G 
is transformed as 

rp (G x \_ ( cos9 sm9\ ( G x \ 
\ G v / V - sind cosd J V G v J ' 

where G x and G y are the G components in the original coordinate system. Therefore under this operation, we have 

T r [S(G)] = S(T r [G]). (12) 



3. Reflection 

The reflection operation is to reflect the system about a line. This operation will not change the form of the 
structure factor. Rather it just change the reciprocal vector. By analogy with the rotation operation, under the 
reflection operation we have 

T re/ [S(G)]=S(T re/ [G]). (13) 



4- Scaling 



The scaling operation includes dilatation and contraction of the coordinates. Consider the operation 
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(14) 



Under this operation 



T s [G-r] = G-T s [r] =T s [G]-r, 



(15) 



and 



dA — > X x X y dA. 



(16) 



From (|T|) and (JTe|) we find 



T S [5(G)] = A !C A 1/ S'(T S [G] 



(17) 



5. The structure factors of circular and elliptic cylinders 

As an example, we calculate the structure factors of a circular and an elliptic cylinder. The latter is considered as 
a result of the scaling operation on the former. For a circular cylinder of radius a, we have 

5(G) = 2TO ^,f a) . (18) 

For the elliptic rod, suppose the lengths of the principal axes (along x and y axes respectively) are a and b. According 
to the idea before, we can obtain the result by a scaling operation on the circular case. Here 

X x = l, >^y = K (19) 



so 



where 



= 2nbMGa) 
AG 



G = dGl + (-fGl (21) 



B. Structure factor of a polygonal cylinder 

Now we consider the general case. Consider a polygonal cylinder with N sides. The N corners of the polygon are 
labeled as rj., rj, . . . , rjv, where the jth corner = (xj, yj), and we define r/y+i — r±. Suppose we can find a vector 
field F such that 

(V x F) z = e~ iGr , (22) 

then according to Stokes' Theorem, we have 

V x F • dA = I F • dv, dA = dAz, (23) 
Jc d 

and thus the original integral (^|) is reduced to computing the line integral 

1=1 F ■ dr, (24) 
Jc d 

where d is the boundary of the polygon. In the rest we present a detailed calculation of this line integration. 
To begin with, assume G x ^ 0. Choose 

F = F y y, (25) 

then solving the equation 

(V X P), = 8 x F y = e- %G r = e -iG x x e -iG y y 

will gives us the solution 

F y = ^e- iG -. (26) 

The integral I now becomes 

• n r ( Xj+1 .y j+1 ) 

I=-^rJ2 e-^ G * x+G yyUy, (27) 

Gx j=1 J{ Xj , yj ) 

where J^-y')^ denotes the line integral between j-th and (j + l)-th corner, i.e., the j-th side. For the j-th side 
we have 

e -t{G x x+G y y)^y 

ijVj+i - yj)(e~ lG r3+1 -e~' G r 
G ' ( r i+i - r i) 

(28) 



Substituting ( |28| ) into (|27| ) and defining 



r i+ i + rj 



(29) 



we find 



A Vj = yj+i-Vj, (3°) 



" zA^e-' G ^sin(G.S J ) 



3=1 



When G K = 0, Eq.(|3l|) no longer holds. We have to find another expression. In this case, since G y ^ 0, we can 
choose 

F = F x x (32) 

and solve 

(V x F) z = -d y F x = ( ' <: ' ( 

to find the solution 

F x = --^e- iG *. (33) 

Taking this F into Eq.(|24|) and defining 

Axj = Xj+i — Xj, (34) 

we find 

/ " -A,y mu:CI-Sj (35) 

. ' Gy G ■ Sj 

3 = 1 J J 

In fact, both ( |3~l] ) and (|35|) can be written in the unified expression 

2iSj -fi 2 e- iG - c ^ sin(G-Sj) 



G • h\ G Sj 

_ A 2^z ■ (ni x S, i, (! <' sin (G ■ Sj) 

where hi is an arbitrary unit vector, and hi is another unit vector defined by hi = z x hi. 
If we choose hi as 



G 

~G' 



where G = |G|, we then obtain a general expression 



3 = 1 



(37) 



^ 2»i-(GxS J )e-' G - c J sin(G.S J ) 

7 -2. G2 erg— ' (38) 



for either G x ^ or G 9 ^ 0. 

In the above we have obtained the general expression Eq. (|3^) of the integral / when G/0. Several remarks are 
worthwhile. 



1. If G = 0, the integral in (||) is simply A, i.e., the cross section area. However, we do not need to consider this 
situation because the structure factor at G = does not involve in the band structure calculation, referring to 
Eq. (|). 



2. Formula ( |31|), (p5|), ( J36| ) and (|3§|) are still correct when the polygon region has several polygons. In this case 
the sum in(p8|) can be rewritten as 

A ^ 2iz ■ (G x S,- )e~ iG - c ^ sin (G ■ S 7 - ) 

P=U P =l Jp 

where p = 1, . . . , P and P is the total number of polygons; j p = 1, . . . , 2V P , with iV p being the total number of 
corners of the p-th polygon. 

3. Since the integral ^3) is calculated along a closed loop, we can always gauge transform the field F by adding 
an arbitrary gradient: 

F' = F + VA, (40) 

where A is a scalar field. Under this gauge transformation, the structure factor is not changed. For example, 
denote F in ( |26| ) as Fi and the F in ( |33| ) as F 2 , then 

/_ e -iG-r\ 

F '- F *= v («d- (41) 

Therefore Fi and F 2 differs by a gradient, and both gave the same structure factor. 



4. If G x = 0, one can use Eq.fl35|); if G y = 0, one can use Eq. (|3l|). It can be proved that the result with G x = 
can also be obtained from Eq. (^Tj) by taking the G x — > limit. Similarly, when G y = 0, one can obtain the 
correct result from Eq. ([H]) by taking the G y — > limit. 

5. For later convenience, we define a function Q(x) as 

. since 

Q{x) = . (42) 



IV. SOME ANALYTIC EXAMPLES 

Now we apply the above formulas to calculate explicitly the structure factors of various configurations. 

A. Square lattice 

First we consider the square lattice. The lattice constant is supposed to be d. The base vectors are 

ai = dx, &2 = dy, (43) 

so the area of a unit cell is 

A = d 2 . (44) 

The bases of the reciprocal lattice are 

bi = -j-x, b 2 = —y, (45) 



which give 



or 



2tt 

G = nibi + n 2 b 2 = — (n 1 x + n 2 y) , (46) 
d 

0.-^0.-^. (47, 



1. Square lattice with circular cylinders (SC) 



First we arrange the circular cylinders in the square lattice. This is simplest type of photonic crystals |29|. The 
radius of the cylinder is a. The structure factor is given by 



5(G) = 2/ 



JxjGa) 
Ga ' 



The filling fraction / and its maximum value f max are 

,2 , „ N j 



/ 



.4 



For a general (raj., n-x) pair we have 



Ga 



\l^-(nl+n 2 2 ) = ^nf(nl+nl). 



(48) 



(49) 



(50) 



2. Square lattice with square cylinders (SS) 

Now we place the square cylinders on the square lattice |?3|. The side length of the square is a. The structure 
factor is given by 



s(G) = fQ{—)Q{-t)=fQ {—) {.—) 



so the area of a unit cell is 



The bases of the reciprocal lattice are 



which give 



(51) 



Here the filling fraction and its maximum value are 

/=(^) 2 , f max = l- (52) 



B. Triangular lattice 

In this case we still assume the lattice constant I and circle radius a. The base vectors are given by 

1 \/3 

ai = dx, a 2 = d(-x + —y), (53) 



A=^-d\ (54) 



i 2tt 1 2vr 2 . 

bl = T (x ~ Ti vl b2 = T7! y ' (55) 



G = mbi + n 2 b 2 = ^ (nix + ^-jJ- -y ) . 



1. Triangular lattice with circular cylinder (TC) 
This is the case considered by [i^j . The radius of the circle is a. The structure factor is 

A Ga 



(58) 



with 



2V3' 



(59) 



For a general (m, 712) pair we have 



'167r 2 a 2 

Ga = \l — ^-^ — (n\ + n\ — nin 2 ) 



3d 2 



(60) 



2. Triangular lattice with hexagonal cylinder (TH) 



In this case we assume the hexagon side length is a, thus the cross section area is 

.2 



A = 



3V3 a" 



(61) 



The filling fraction is given by 



A /a\ 2 



(62) 



with maximum value 



The structure factor is 



(63) 



5(G) = MS Ln 
G x a 



(3G X + V3G y 



(G x — \/3Gy 



sin 



(3G K - V3G y )a 



Q 



(G x + ViG y )a 



(64) 



S (G) = M^ M11 



niira 



(ni + n 2 )ira 



(rii — ri2)iTa 



sin 



(2ni — 7^2)™ 



ri27ra 



(65) 



C. Tables of structure factors 



In this section we consider various arrangements of rod cross section and lattices. Eighteen patterns (configurations) 
are illustrated in Fig. ([!]) and (J2J) . Here we consider square, triangular, and honeycomb lattices. The abbreviation of 
these patterns are defined in the figure captions. For example, TC refers to the situation that circular cylinders are 
arranged in a triangular lattice. For a circular rod the parameter a denotes the radius of cylinder. For a polygonal 
cylinder, a represents the side length. For a "cross" -shaped or a "rotated cross" -shaped rod, a and era represent 
the "long" and "short" side lengths of the corresponding rectangle. In a unit cell of the "triangular-diamonds" or 
"triangular-rotated diamonds" pattern (See Figs. 0(j) and (k)) there are three "diamonds" . In the "rotated diamonds" 



pattern the distance between the center of each diamond and the center of unit cell is chosen to be d/4; d is the lattice 
constant. 

The crystal structures shown by these two figures are most common in photonic design. While some have already 
been published, many are reported for the first time. In the next section, a new class of photonic crystals will be 
investigated. The corresponding band structure will be computed. The corresponding properties and the forms of 
unit cells are given in Table I, and the structure factors are listed in Table II, respectively. In these tables, the first 
column refers to the combination of the lattice structure and the shape of rods, the second to the unit cells, the third 
to the filling factor formula, with the maximum value being listed in the fourth column. The next two columns refer 
to the two base reciprocal vectors. 



V. THE NUMERICAL RESULTS 



In this section we consider the band structure of the 2D photonic crystals made by drilling hollow cross-shaped 
cylinders in an alumina ceramic background and by placing cross-shaped alumina ceramic cylinders in the air. The 
cross-shaped cylinders are arranged in a square lattice configuration, corresponding to Fig. |l](d). The following 
parameters are used: filling factor / = 0.5, the ratio between b and a is 0.3. The dielectric constant of the alumina 
ceramic is 8.9. We used more than 600 plane waves. The inaccuracy is less than 1%. 

The photonic band structures are plotted in Fig. ||. Here we see from (a) that an absolute band gap appears in the 
high frequency range for the E- and H- polarizations, in the case of air cylinders in the dielectric medium. Incidently, 
the absolute band gap position is at around the same regime as in |Q and |Q ; we note that in (iH) , it was the very 
high dielectric cylinders that are placed in air. This indicates that it is relatively easy to obtain the absolute band 
gap in the infrared regime with this type of photonic crystals. We also see that there are more complete gaps in the 
E-polarization mode than in the H-mode. These occur at around u>d/(2irc) = 0.25,0.6, but the largest gap happens 
at the position of the absolute gap around u>d/ '(2irc) — 0.7, where the gap of the H-polarization is wider than that of 
the E-polarization. In the opposite case where the cylinders are made of high dielectric constant materials, there is 
no absolute band gap at all. In fact, there is no gap for the E-polarization. However, a complete gap appears around 
ujdj (2irc) — 0.35 for the H-polarization. 



VI. SUMMARY 



In this article, we derived necessary formulas needed for computing the band structures of various types of photonic 
crystals. The results are systematically presented in tables, so that the reader could easily make a use of them. As an 
example, we applied the results to two dimensional arrays of hollow cross-shaped cylinders embedded in an alumina 
ceramic medium. For this new class of photonic crystals, a large absolute phontonic band gap is discovered in the 
high frequency region. 

Finally, we want to make a connection between the photonic and sonic crystals. In the acoustic case, the wave 
equation is governed by 

V ' (W) VP(r) ) + iSiif^ = ^ (66) 
In the H-polarization mode, Maxwell equation is from Eq. (|l|) 

V '(^ii(')) + ^)^W = °. ^ 

where c 2 = Cg/e(r). Therefore the acoustic pressure field and the magnetic field have an one-to-one correspondence 
given by the following mapping, 

p(r)^e(r), p(r)c 2 (r) ^ e(r)c 2 (r). (68) 
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FIG. 1. Patterns of the square lattice, (a) square-circle (SC). (b) square-square (SS). (c) square-rotated square (SRS). (d) 
square-cross (SCR), (e) square-rotated cross (SRCR). (f) square-mixed square (SMS), (g) square- mixed cross (SMCR). (h) The 
first Brillouin zone of square lattice. The band structure is calculated along the boundary M — T — X — Mof the gray region. 
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FIG. 2. Patterns of the triangular lattice, (a) triangular-circle (TC). (b) triangular- hexagon (TH). (c) triangular-rotated 
hexagon (TRH). (d) triangular-triangle (TT). (e) honeycomb-circle (HC). (f) honeycomb- hexagon (HH). (g) honeycomb-rotated 
hexagon (HRH). (h) honeycomb-triangle (HT). (i) honeycomb-rotated triangle (HRT). (j) triangular-diamonds (TD). (k) tri- 
angular-rotated diamonds (TRD). (1) The first Brillouin zone of triangular lattice. The band structure is calculated along the 
boundary X— T — J — Xof the gray region. 




FIG. 3. The photonic band structures. The solid and dashed lines refer to the E-polarization and H-polarization modes 
respectively, (a) Air cylinders in alumina ceramics. The shaded area denotes the complete band gap. (b) Alumina ceramics 
cylinders in air. 
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TABLE II. Structure Factors for Square and Triangular Lattices 
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